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In quantum physics, measurements give random results and yield a corresponding random back 
action on the state of the system subject to measurement. If a quantum system is probed con¬ 
tinuously over time, its state evolves along a stochastic quantum trajectory. To investigate the 
characteristic properties of such dynamics, we perform weak continuous measurements on a super¬ 
conducting qubit that is driven to undergo Rabi oscillations. From the data we observe a number 
of striking temporal correlations within the time dependent signals and the quantum trajectories of 
the qubit, and we discuss their explanation in terms of quantum measurement and photodetection 
theory. 


A collection of quantum oscillators prepared in the 
same initial state will oscillate in phase until decoherence 
mechanisms cause the ensemble to dephase. If the oscil¬ 
lators are weakly measured, the resulting average signal 
will reflect the ensemble behavior of the oscillators and 
exhibit damped oscillatory behavior. In addition, if the 
quantum oscillator is not in an eigenstate of the mea¬ 
sured observable, additional dephasing will be present as 
mandated by the Heisenberg uncertainty principle. Thus, 
stronger measurements result in faster damping of the ob¬ 
served oscillatory signal. But, what average signal does 
one expect to observe if instead a projective measurement 
is used to post-select [T] the oscillators in the same final 
state? In this Letter, we show with experimental data 
that the average post-selected signal is exactly the time 
reverse of the pre-selected one where the oscillators are 
prepared in the same initial state. This temporal cor¬ 
relation between qubit state populations and observed 
measurement signals is an example of the interwoven na¬ 
ture of measurement signals and quantum trajectories 
that have become available with current experiments on 
superconducting qubits W- Inspired by their use in 
the analysis of stochastic processes, we address different 
two-time correlation functions, available from the exper¬ 
iments and discuss how their properties may be analyzed 
by quantum measurement theory and by master equation 
approaches. 

In quantum measurements, the detection of light con¬ 
veys information about the state of the emitter and in¬ 
duces corresponding backaction on its quantum state. 
Temporal correlations of light reveal purely quantum ef¬ 
fects and have been essential to the development of quan¬ 
tum optics 00 which traditionally applies the master 
equation and the quantum regression theorem for eval¬ 
uating the observed measurement signal. Recent exper¬ 
iments have harnessed the quantum trajectory descrip¬ 
tion nmE] to track the emitter’s dynamical evolution. 
In this work, we take such analysis one step further 
and investigate two-time correlation functions between 


the measurement signal and the quantum state. As ex¬ 
pected from quantum optics, the correlations associated 
with quantum measurements at different times can re¬ 
veal purely quantum effects for example by testing the 
assumptions of macrorealism via the Leggett-Garg in¬ 
equality mm- 

We perform continuous weak measurement on a super¬ 
conducting qubit in the energy basis while the qubit is 
resonantly driven to produce Rabi oscillations go snug. 
Measurements are executed in the dispersive regime of 
cavity quantum electrodynamics m where the inter¬ 
action between the qubit and a mode of the cavity is 
given by the interaction Hamiltonian, H lni = —HxcP cl<t z , 
where a'(a) is the creation (annihilation) operator for the 
cavity mode, a z is the Pauli pseudo-spin operator that 
acts on the qubit in the energy basis and x is the dis¬ 
persive coupling rate. Homodyne probing of the cavity 
can thus be used to conduct both weak and strong mea¬ 
surements of the qubit state in the a z basis. We describe 
these measurements by the theory of POVMs (positive 
operator valued measures) which relates the homodyne 
voltage signal V to the qubit state by the operators im 

Ul v = (27ra 2 )- 1/4 e- (v -' T d 2 / 4 “ 2 . (1) 

The probability of detecting a homodyne voltage V, at 
time t, P(Vt) = Tr(fiyptf2y), depends on the density 
matrix p t , and it is the sum of two Gaussians centered 
at V = ±1 and weighted by the qubit state populations 
Pt° and p\ l . Depending on the variance a 2 = \/{Akp/Si), 
where k represents the measurement strength, p is the 
quantum measurement efficiency and At is the integra¬ 
tion time of the measurement, this operator describes 
both strong (projective) and weak measurements of the 
qubit. 

If the measurement is strong, the variance a 2 is small 
and the measurement outcome unambiguously belongs 
to one of the two disjoint Gaussian distributions [T8j - If, 
conversely, the measurement is weak, the variance a 2 is 


2 


large, and P(V) can be approximated by a single Gaus¬ 
sian distribution centered at the expectation value of er z , 

P(V) ~ (2tt a 2 )- 1/2 e- ( - v -^ ))2/2a2 . (2) 

In this regime, V thus provides a noisy estimate of (a z ). 

In addition to the assignment of outcome probabilities, 
the POVM operators describe the effect of the measure¬ 
ment on the quantum state, p —> ClypSVy. In the limit 
of strong measurements, a result V ~ ±1 corresponds 
to the projection operators 0± 2 = [cr z ± l)/2 on the ±z 
qubit eigenstates. For weak measurement, the backaction 
on the state is small and the evolution of the density ma¬ 
trix associated with a Hamiltonian Hu (describing uni¬ 
tary evolution) and a measurement signal V(t) can be 
obtained from the stochastic master equation in mum, 

P = [Hr, p] + k( a zPVz - p) 

+2rjk(a z p + pa z - 2Tr(a z p)p)V(t). (3) 

If we disregard the measurement outcome, or average 
over its values, the qubit density matrix is still subject to 
dephasing due to the probing and solves a deterministic 
master equation, 

p\det = ~p[H Rl p} + k{a z pa z - p). (4) 

in 

In our experiment, the qubit is formed from the two 
lowest energy levels of a transmon circuit [2D] (see Fig. 
1(a)). The transmon is clispersively coupled to a three 
dimensional waveguide cavity |21] . The qubit transition 
frequency is uj q /2TT = 4.0033 GHz and the cavity tran¬ 
sition frequency is 6.9914 GHz with bandwidth k/2i r = 
9.88 MHz and dispersive coupling rate x/27t = —0.425 
MHz. Since 2|\| <C k, qubit state information is encoded 
in a single quadrature of a transmitted microwave tone 
that probes the cavity. Qubit decay from the excited to 
the ground state is slow and can be ignored for the time 
scales of our experiments. The reflected signal is am¬ 
plified by a near-quantum-limited Josephson parametric 
amplifier [122] and demodulated and digitized after fur¬ 
ther amplification at room temperature. Further details 
of the experimental setup can be found in supplemental 
information [23] . 

We prepare the qubit in an initial state by applying 
a 7r/2-rotation about the y axis and then making a pro¬ 
jective measurement This projective measurement 

can be used to herald either +z or — z as the initial state 
or, if the result of the projective measurement is ignored, 
to prepare an initial mixed state. Following this prepara¬ 
tion, the qubit is subject to continuous rotation given by 
Hr = h£ln<j y /2, where Qr/2tt = 1.16 MHz is the Rabi 
frequency, and probing as given by the operators fly. 
The experiment is repeated several times to form an en¬ 
semble of measurement signals {Vi(t)} which is analyzed 
in a post processing step. 



Figure 1: (a) Experimental setup consisting of a transmon 

circuit coupled to a 3D cavity, (b) Continuous probing of 
a qubit driven with a Rabi frequency Q.u/2-k = 1.16 MHz 
after its heralded preparation in +z or —z. The probe is 
given by the measurement operator \/ka z , where k = k, 

parametrizes the measurement strength (k/2iv = 95 kHz) and 
rj = 0.35 is the quantum measurement efficiency, (c) Pre¬ 
selected average of the output measurement signals, condi¬ 
tioned on the heralded preparation of +z. (d) The expecta¬ 
tion value Zi(t) in a number of individual quantum trajecto¬ 
ries (gray lines) and their mean value (blue line), which is in 
agreement with the mean signal shown in the upper panel. 

Fig. 1(b) displays how we obtain the measurement sig¬ 
nals and how the qubit is initially heralded in the +z 
state by a strong measurement. Using the stochastic 
master equation © we propagate the density matrix 
according to the individual, weak measurement signals. 
The average voltage signal V ( t ) in Fig. 1(c) shows damped 
Rabi oscillations corresponding to the gradual dephasing 
of the qubit due to the measurement interaction. In Fig. 
1(d), the qubit expectation value, Zi(t) = Ti'(pi(t)a z ) cal¬ 
culated by solution of © , is shown for several of the indi¬ 
vidual quantum trajectories HUM], confirming that the 
damping of the measurement signal is due to dephasing 
of the ensemble. 

In Fig. 2(a), the black curve displays the average of the 
measurement signals conditioned on the outcome of a fi¬ 
nal projective measurement which post-selects the qubit 
in the +z state. The resulting average exhibits an oscil¬ 
latory signal that is damped backwards in time [4[. Sur¬ 
prisingly, this is exactly the time reverse of the damped 
signal observed in Fig. 1(c) and it exhibits the same full 
contrast at the final time t = T that we observe for the 
pre-selected average at t = 0. 

Fig. 2(b) displays a sample of the trajectories z,;(f), 
which are propagated forward from the initial state using 
the stochastic master equation. Since the post-selection 
average is conditioned only on the final measurement, a 
roughly equal number of the trajectories, included in this 
average, are initially detected in the +z and —z states. 
As shown in Fig. 2(c) immediately prior to the post¬ 
selection, the trajectories take on several different val- 
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Figure 2: (a) The average of measurement signals conditioned 
the post-selection of +z (black trace). The red trace is the 
average of all of the measurement signals weighted by the fi¬ 
nal population p°°(T). (b) The quantum trajectories whose 
measurement outcomes contribute to the post-selected and 
weighted averages in (a) are shown as Zi(t) evolves forward 
in time from the initial state, (c) The magnified view of 
the last time segment shows that the trajectories populate 
the qubit eigenstates very differently. Trajectories that are 
post-selected in the +z state are indicated with the open cir¬ 
cles. The measurement signals Vi corresponding to this sub¬ 
ensemble form the average Fps- (d) The vertical bars indicate 
the relative weighting applied to all of the measurement sig¬ 
nals that contribute to Fwp- 


ues of Zi(T ) between ±1 , and thus it may be surprising 
that the average measurement signal (black curve in Fig. 
2(a)) exhibits full contrast oscillations at the end of the 
sequence. 

It is tempting to assume that the time reversed damped 
signal is merely a quirky consequence of post-selection 
[TJ [25], but we will now show that the same behavior can 
be obtained by a suitably weighted average over all tra¬ 
jectories without post-selection. Consider the stochastic 
master equation © yielding for each experimental run 
an independent realization pi(t) of the time dependent 
conditioned density matrix of the qubit. We now aver¬ 
age all the measurement signals Vi(t), weighted according 
to the final state probability P(i, +z ) = p°°(T), 


V wp 


E iPTirmt) 
EiP°°( T ) 


(5) 


Fig. 2(d) indicates using vertical bars the relative contri¬ 
bution of each of the signals to the average. As shown 
by the red curve in Fig. 2(a), this weighted average is in 
excellent agreement with the post-selected average. 

In order to explain the equivalence between the post- 
selected and the weighted average we observe that the 
post-selected mean value can be written as, 


l^PS 


E i«i 


( 6 ) 


where the stochastic variable rii = 0,1 depends on the 
post-selection result for the i th element of the ensem¬ 
ble. Since the probability of the post-selection mea¬ 
surement yielding +z is p°°(T), in the limit of a large 
ensemble, Ei n *E(t) approaches JT/9°°(T)V)(f). Simi¬ 
larly, the denominator n,; approaches V ; p°°(T), and 
Vwp = Vps because the individual trajectories faith¬ 
fully predict the probability of the final projective state 
measurement. 

To explain the symmetry of the pre- and post-selected 
averages, let us consider the joint probability P(Vt , ±zt) 
for the measurement outcome V) at time t and a pro¬ 
jective measurement on ±zt at time T. This probability 
factors as the product of the probability for the first mea¬ 
surement outcome P(V t ) and the conditional probability 
P(±ZT\V t ), 

P(V u ±z T ) = P(V t )- Tr{ e ^-*)[^^]| ± z)(±z|},(7) 

where p is the state prior to the weak measurement of V), 
Vlyptiy/P(Vt) yields the normalized state conditioned on 
the outcome Vt, and e L ( T ~ t ' ) denotes the linear propaga¬ 
tor from time t to T of the density matrix according to 
the master equation ©• 

It is convenient to express the deterministic mas¬ 
ter equation evolution from time t to T as a Kraus 
map, e L ^ T ~ tS) \p] = E« K a pK} y , with operators obeying 
e a KK a = i mm- By making use of the cyclic 
properties of the trace, we can shift the Kraus operators 
K a to the right hand side of the expression |7]), yielding, 

P(V t , ±z T ) = Tr {n v p(t)n{. J2 K l\± z){±z\K a } 

a 

= Tr [n v p(t)nlE(t)}, (8) 

where E(t) is defined as the operator expression involving 
the final state | ± z)(±z| and the adjoint Kraus map. We 
do not need the explicit form of the Kraus operators, as 
E(t) can be found by solving the (adjoint) master equa¬ 
tion, 

\W,>,E\+t^Ea 2 -E), ( 9 ) 

backwards in time from the state E(T) = | ± z)(±z| to t. 

While it may seem unusual to discuss the probability 
and weighted average of a measurement outcome V) con¬ 
ditioned on a later projective measurement, such anal¬ 
yses have recently been proposed m and successfully 
applied to experiments IB SI ESI. If we do not condition 
on the final measurement and set E(t) oc I in Eq.®, we 
recover the usual prediction, V = (a z ) = p 00 (t ) — ? 1 (t), 
while if we condition on the outcome of the last and not 
the initial a z measurement, p is proportional to the iden¬ 
tity matrix in and Kps = E 00 (t) — If 11 (t). The only 
difference between the forward evolution Q of pit) and 
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Figure 3: Two-time correlation function between the mea¬ 
surement signal Vi(t) and the inferred qubit density matrix 
element p°°(t'). The side panels show the correlations as func¬ 
tion of t for fixed values of t!. 


the backward evolution ^ of E(t) is the sense of rotation 
of the damped Rabi oscillations. They therefore yield 
identical heralded predictions (see analytical expressions 
in the supplementary information) for the voltage signal. 

While mysterious action from the future through post¬ 
selection evokes fascinating scientific debate [I] [25 ] [29l [30] 
the predictions we make for the post-selected or weighted 
averages merely reflect the correlation between the qubit 
observables at different times. Such correlations are 
central in the quantum optical characterization of light 
sources, and the time evolution between times t and T in 
our Eq.([8]), indeed, appears in a very similar manner in 
the quantum regression theorem m when it is applied to 
calculate field intensity and amplitude correlations [32|- 

To examine further the relationship and the tempo¬ 
ral correlations between the measurement signal and the 
inferred qubit density matrix, we show in Fig. 3 the "hy¬ 
brid" correlation function between the measured signal 
and the inferred qubit state, 


Vwpitjt') 




( 10 ) 


obtained as an average over all the experimental runs of 
the voltage signal Vj(t) at time t weighed by the qubit 
excited state population Pi°(t') at time t! . 

The side panels in figure 3 show V\yp{t,t') as a func¬ 
tion of t between 0 and T for fixed values of t!. We clearly 
observe a change in the behavior of the function at t = t! 
which suggests a different correlation regime before and 
after t'. This has a natural explanation: V)(f) is the sum 
of a term proportional to ( a z ) and a white noise term 
W(t). When t < t ', the noise term W(t) contributes a 
change of the conditioned density matrix due to ^ hence 


it affects all later values of p^(t'). We thus expect that 
until t = t' , the product in Eq. (101 will contain a con¬ 
tribution that is quadratic in the white noise term W(t ), 
while for t > t', W(t) is uncorrelated with the earlier 
value and hence it averages to zero in the sum 

over i in Eq. (|l0| ). 

For t < t' , Vwp{t,t') can be calculated by the analy¬ 
sis that we applied to the post selected data (by setting 
T = t') and it exhibits all the features we have explored 
so far. Even though we do not perform a projective mea¬ 
surement at time t', due to the equivalence of Rps and 
R WP , the function plotted in the region with t < t' can 
be calculated from ([8]), with E(t) propagated backwards 
from the state | + z)(+z\ at time t' [23]. 

We recall that p®°(t') = ((cr-(t'))i + 1)/2, and V{t') = 
{a z (i')) + W(t '), where W(t') has zero mean and is uncor¬ 
related with all previous quantities. For t < t! , this leads 
to an alternative expression for Vwp{t-i t') involving only 
the mean value and the two-time correlation function of 
measured quantities, 


Rwp (t, t') = 


ViiWW + Viit) 


V(t>) 


1 


( 11 ) 


Note that the signal-state and the signal-signal correla¬ 
tions being averaged in the different expressions for Fwp 
are not identical and equality between them hold only for 
large ensembles of measurement records. 

It is interesting to observe that for t < t' where the 
state-signal correlation function can be expressed solely 


in terms of experimental signal correlations (111, the past 


quantum state formalism provides a deterministic theory 


for the correlation function (10). This is, indeed, inti¬ 


mately related to the conventional use of the quantum 
regression for the same purpose in quantum optics |32j . 
For t > t ', the correlation function may instead be writ¬ 
ten Vwp{t,t') = Pi -£ Such correlations 

and higher moments of the quantum trajectory solutions 
do not obey any simple deterministic equation, and for 
each possible density matrix p(t') resulting from the pre¬ 
ceding stochastic measurement record ([3]) , we must prop¬ 
agate p{t) by Eq. (|3|) or Q to determine p 00 (t , )p 00 (t). 

Two-time correlation functions like p 00 (t , )p 00 (t) have 
so far not found applications in quantum optics. The 
recent progress in circuit QED, however, has made it 
possible to assess the dynamics of single quantum sys¬ 
tems and thus demands a more detailed characterization 
of the stochastic system dynamics during measurements. 
Recently, a theory was presented for the most likely path 
of a stochastically evolving density matrix El m 133, and 
correlation functions, indeed, present a natural quantita¬ 
tive characterization of such a theory. Finally, stochas¬ 
tic state dynamics is an unavoidable, and even useful, 
component of precision metrology and parameter esti¬ 
mation by continuous probing |35], where the so-called 
Fisher Information is, in fact, given by an ensemble aver- 
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age of products of elements of density matrices that solve 
stochastic master equations |.'lfi| . 
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Figure 4: Experimental setup. The weak measurement is always on except when the projective measurements are performed. 


I. EXPERIMENTAL METHODS AND MEASUREMENT CALIBRATION 


Figure 4 shows a detailed diagram of the microwave setup for the experiment which uses the same qubit and 
amplifier circuits as described in previous work [T]. 

To calibrate the parameters of the weak continuous measurement we perform separate calibration experiments. We 
prepare the qubit in the +z or —z states using a projective herald measurement and then digitize the ensuing weak 
measurement. These weak measurement results are Gaussian distributed and we scale the recorded signal such that the 
distributions are centered at +1 and —1 for the +z and — z states respectively, P(V t \ +z) = (27ra 2 ) -1 / 2 exp[— ^2 ]> 

P(Vt\ — z) = (27ra 2 )~ 1 / 2 exp[— ]• The quantum measurement efficiency 77 = r] co \ij amp depends on both the 

collection efficiency and added noise from the amplifiers and is determined from the variance, a 2 = l/4fc?yAf = 
k/{\ d>x 1 m)/S.t) where k = 4y 2 n/K represents the measurement strength and At is the integration time. The parameter 
k is also related to the characteristic measurement time r [2] as k = 1/4?7t. We used Ramsey measurements to calibrate 
both the ac Stark shift (2 \n) and the dephasing rate (T m = 8x 2 n/n) for photon numbers ranging between n = 0 and 
n = 1.3, thus determining n and %. 
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II. STATE PROPAGATION 

We propagate p forward from the initial state p 00 = 0.95, p 01 = 0 when the initial qubit state is +z and p 00 = 0.05, 
p 01 = 0 when the initial state is — 2 . The stochastic master equation displayed in the main text omits qubit state 
dephasing for clarity, but the quantum trajectories are calculated including this dephasing which is characterized by 
the rate 7 = l/TJ, where TJ = 16 ps. To account for this, k should be replaced by k + 7/2 in the second term in Eq. 
(3), and in Eq. (4) and (9) in the main text. The full stochastic master equation is given by, 

^ = + ( k+ \^ a -P a ~ ~ p) + 2i 1 k(a z p + pa z - 2Tr(cr-p)p)Vt. (12) 

The Rabi frequency is 0^/27r = 1.16 MHz and we use time steps of At = 20 ns. We perform quantum state tomography 
on an ensemble of experimental iterations with similar measurement values to verify that we have accurately inferred 
the quantum trajectory as we have done in previous work m- 


III. ANALYTICAL EXPRESSIONS FOR V W p 


From Eq. ( 8 ) in the main text we can calculate the conditioned expectation value of the voltage signal at time t [T] . 
The joint probability of measurement outcome V) and subsequent measurement outcomes and evolution accumulated 
in the matrix E can be evaluated, 


P p (V t ) oc p00£00 e -W t -l) 2 /2 a 2 + p ll E ll e - {Vt + l) V2a 2 + (p 10^01 + p 01 E 10y- (V?+l)/2a 2 _ 


(13) 


Where the subscript p denotes “past" in that the probability depends on both earlier and later probing results. Thus 
the mean value is given by V p = f P p (V)VdV and can be evaluated. In the limit of weak measurement (a —> 00 ) we 
have, 


y PS = 


p oo E oo _ p u E 


11 


p oo E oo + p U E ll + p 01 E 10 + p io E oi 


(14) 


where the elements of p and E are all evaluated at time t. 

The deterministic master equations for p and E (Eq. (4) and (9) in the main text) can be solved analytically. The 
two components of p are given by 


P°°(*) = e fct (^)[sin(Tt) + ^ cos (Tt)\ + 


(15) 


p 01 (t) = e- kt (^) S in(rt), (16) 

where T 2 = — k 2 and assuming the qubit is in the ground state at t = 0. The analytical solution for the backward 

evolution, Eq. (9), can be obtained as well for the final state +z at t = T. 

E°\t) = e- fe (T-*) ( l. )[sin(r(r _*)) + £ cos(T(T - i))] + (17) 

E 01 (t) = -e- fe ( T -‘)(^) sin(T(T - i)). (18) 


IV. COMPARISON OF PRE-SELECTED AND POST-SELECTED MEASUREMENT SIGNALS. 

Figure 5 shows the averaged measurement signals conditioned on the herald measurement and the final projective 
measurement. Aside from a slight difference at the beginning (end) which arises from transients in the measurement 
strength resulting from the herald measurement, the two curves are in close agreement. 
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Figure 5: Pre-selection and Post-selection average measurement signals. The black curve (top axis) corresponds to the average 
measurement signals from the qubit post-selected in the ground state, which is plotted in reverse. The gray curve (bottom 
axis) is the average measurement signal corresponding to pre-selection in the ground state. 
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